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SUMMARY 


The  work  reported  here  draws  attention  to  the  geometrical  disposition  of 
shield  elements  as  an  Important  means  of  enhancing  shield  efficiency  in  a 
situation  where,  at  first  sight,  the  preponderance  of  Compton  scattering  seems 
to  preclude  improvements  by  any  method  other  than  increasing  the  total  mass. 
The  principal  effort  in  Phase  I  has  been  to  demonstrate  the  practical  utility 
of  using  good  geometry,  for  which  purpose  it  has  been  necessary  to  develop  a 
mathematical  description  of  what  good  geometry  -  hitherto  a  vague  and  qualita¬ 
tive  concept  -  actually  is.  The  advantage  of  layering  the  shield  or  of  making 
other  geometrical  modifications  relying  on  the  same  theme,  such  as  introducing 
thickness  gradients,  has  also  been  investigated  for  situations  in  which  the 
theoretically  ideal  geometry  cannot  be  achieved  because  of  fabrication 
constraints. 

The  fact  that  dose  buildup  factors  can  be  as  high  as  10  or  even  100 
indicates  the  tremendous  contribution  that  scattered  radiation  can  make  to  the 
overall  dose  on  the  far  (safe)  side  of  a  gamma-ray  shield.  The  layering 
principle  seeks  to  reduce  this  contribution  by  configuring  the  shield  in  such 
a  way  as  to  allow  increased  leakage  of  the  scattered  photons  from  the  shield. 
In  the  course  of  studying  this  principle,  a  powerful  and  versatile  Monte  Carlo 
formalism  has  been  developed  which  is  easily  extended  to  allow  investigation 
of  other  more  complex  shield  configurations  involving  for  instance  density 
gradients,  lenticular  layer  profiles  and  curved-  or  conic-shaped  shield 
members. 

The  Importance  of  this  Phase  I  study  Is  that  It  verifies  the  dramatic 
effect  of  shield  configuration  on  target  dose,  indicating  the  potentiality  of 
methods  for  shield  design  other  than  the  conventional  mass-increase  approach, 
and  thus  illustrates  the  need  to  include  geometric  optimization  as  an  essen¬ 
tial  component  In  overall  shield  design.  In  this  respect,  a  variant  of  the 
Monte  Carlo  method,  called  Inverse  Monte  Carlo,  may  be  useful  In  future 
studies. 
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I.  INTRODUCTION 


Reduction  of  shield  mass  is  of  the  greatest  importance  in  the  shielding 
of  both  space  and  terrestrial  vehicles  -  in  the  space  situation  in  order  to 
minimize  lift-off  weight  and  in  the  terrestrial  case  because  of  terrain- 
loading  and  motive-power  constraints.  Efforts  to  reduce  the  mass  are  severe^ 
complicated  in  the  range  of  gamma-ray  energies  delivered  by  nuclear  weapons  by 
the  fact  that  Compton  scattering  is  the  principal  attenuation  mechanism  in 
that  energy-range.  In  these  circumstances  the  attenuation  in  a  shield  of 
given  cross-sectional  area  is  independent  of  the  nature  of  the  absorbing 
material  (with  the  exception  of  hydrogen)  and  thus  the  total  mass  is  the 
primary  factor  determining  shield  efficiency. 

This  preliminary  study  examines  the  possibility  of  using  the  properties 
of  the  Compton  scattering  process  to  Improve  shield  performance  for  a  given 
mass,  or  to  reduce  the  mass  required  for  a  specified  performance.  The  two 
results  of  a  Compton  collision  are  (1)  partial  energy  loss  ("degradation")  and 
(2)  deflection  ("scattering")  away  from  the  precollision  photon  trajectory 
which,  for  first  collisions,  is  the  llne-of-sight  connecting  the  source  and 
target.  It  is  clear  that  for  a  sufficiently  large  distance  d  between  the 
initial  collision  site  and  the  shielded  object  ("target"  in  the  sequel)  an 
arbitrarily  large  radial  displacement  r  from  the  1 ine-of-sight  can  be  achieved 
if  the  photon  collides  only  once  in  traversing  the  shield  (Fig.  1-1).  Indeed, 
substantially  all  of  the  scattered  photons  emerging  from  a  shield  can  be  made 
to  miss  the  target  provided  the  ratio  d/R  (where  R  is  shield  radius)  is 
sufficiently  large  and/or  the  shield  is  thin  enough  so  that  multiple  scat¬ 
tering  does  not  affect  the  attenuation  -  if  there  is  more  than  one  collision, 
a  scattered  photon  that  would  otherwise  have  missed  the  target  may  be  deflect¬ 
ed  inwards  towards  the  1 ine-of-sight  once  again. 

The  Fig.  1-1  situation  is  usually  referred  to  as  "good"  geometry  while 
the  case  of  very  thick  shields  and/or  small  d/R  ratios  (Fig.  1-2),  where  the 
potential  shielding  ability  of  the  absorber  is  reduced  by  multiple  scattering, 
is  called  "poor"  geometry  [1].  In  this  Phase  I  study  we  have  attempted  (1)  to 
quantitatively  investigate  the  notion  of  "good"  geometry,  and  (2)  to  test  a 
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novel  modification  of  it  in  which  splitting  or  layering  of  a  shield  mass  may 
improve  performance  by  a  repetitive  use  of  good  geometry.  The  application  of 
these  and  related  concepts  to  practical  situations  is  beyond  the  scope  of  this 
investigation  but  will  be  addressed  in  a  Phase  II  submittal. 

A.  Good  Geometry 

As  far  as  the  Investigators  are  aware,  no  deliberate  attempt  to  make  use 
of  good  geometry  (much  less  of  the  modification  of  it  we  propose)  to  enhance 
shielding  has  been  made,  although  use  of  pod  geometries  is  sometimes  forced 
upon  one  in  attempts  to  achieve  sore  unrelated  goal,  e.g.,  to  give  consistent 
results  in  certain  absorption  meastrements  [2]  or  as  a  means  of  improving  film 
co'  trast  in  radiographic  magnification  studies  where  the  large  d/R  ratio  is 
usod  to  reduce  the  effects  of  scattered  radiation  on  the  sensor  [3];  there, 
however,  the  total  dose  to  the  sen  or  is  essentially  unaltered  and  the  dose  tc 
the  target  (i.e.,  the  patient)  is  m  fact  increased. 

The  calculational  method  used  in  dealing  with  radiographic  magnification 
is  semianalytic  and  of  very  limitec  power.  The  alternative  Monte  Carlo  (MC) 
method  is  much  more  versatile  and  powerful  by  virtue  of  its  ability  to  treat 
situations  of  almost  arbitrary  complexity.  It  has  been  applied,  for  example, 
to  the  calculation  of  dose  or  flueree-related  quantities  at  or  immediately 
contiguous  to  the  exit  face  of  plane-parallel  slabs,  such  as  in  the  determina¬ 
tion  of  build-up  factors  [4],  but  does  not  seem  to  have  been  applied  to  the 
calculation  of  dose  quantities  over  the  distant  planes  in  which  good  geometry 
Is  achieved.  It  was  therefore  proposed  to  develop  the  MC  formalism  In  order 
to  study  good-geometry  conditions  and  the  layering  principle. 

B.  The  Layering  Principle 

The  premise  of  the  layering  proposal  is  that  the  exit  face  of  a  thick 
shield  is  "fed"  by  a  photon  population  comprised  of  (1)  a  relatively  small 
number  of  primary  (undeflected)  photons  accompanied  by  (2)  a  much  larger 
number  of  scattered  photons.  Just  as  the  arrival  of  the  scattered  exit  flux 
at  the  target  can  be  considerably  reduced  by  using  good  geometry,  the  "feed" 
of  multiple-scattered  photons  to  the  scattering  layer  immediately  adjoining 
the  exit  face  can  be  reduced  by  ensuring  that  the  last  layer  Is  in  good 
geometry  with  respect  to  the  preceding  layers,  i.e.,  by  layering  the  absorber 
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(Fig.  1-3).  8y  introducing  one  or  a  series  of  similar  layers,  the  feed  from 
any  one  layer  of  the  absorber  to  its  successor  can  be  reduced;  thought  of  in 
another  way,  a  layered  geometry  gives  the  photons  scattering  in  a  given  layer 
an  opportunity  to  "leak"  laterally  (and  thus  escape  collision  in  deeper  layers 
or  in  the  target),  while  in  an  unlayered  shield  of  the  same  total  thickness 
lateral  leakage  is  partially  compensated  for  by  scattering  toward  the  line-of- 
sight.  The  thrust  of  the  Phase  I  proposal  to  investigate  these  matters  was 
therefore  to: 

(1)  Develop  the  MC  formalism,  which  involves  the  preparation  of 
cross-section  libraries,  transport  models  for  the  specific 
geometry  considered  and  techniques  for  scoring  the  desired 
quantities.  This  permits  the  conditions  for  the  existence  of 
good  geometry  to  be  ascertained  for  any  photon  energy  and  any 
shield  material . 

(2)  Set  up  a  code  capable  of  dealing  with  arbitrary  layer  thick¬ 
nesses  and  separations,  shield  materials  and  input  spectra,  and 
verify  the  operation  of  the  code. 

(3)  Use  the  code  to  test  the  improvement  In  the  effectiveness  of  a 
given  shield  mass  brought  about  by  adjusting  its  configuration, 
i.e.,  to  study  the  geometrical  effects  of  shielding. 

The  initial  code-implementation  and  application  phase  has  been 
successfully  completed  and  is  described  in  this  report. 
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(a)  Contiguous  Layers 


(b)  Separated  Layers 


A  photon  scattered  in  Layer  1  can  be  redirected  toward  the  target  if  it 
enters  Layer  2  (a),  but  "leaks"  from  the  system  when  the  layers  are 
separated  (b). 


Figure  1-3.  The  Effects  of  Layering 


II.  application  and  description  of  the  layered  shield  geometi. 


A .  Potential  Applied  :> is 

The  concept  of  gam  .1  -ray  shield  layering  and  configuration  design  has 
potential  military  applications  for  the  obvious  reasons  that  reduced  shield 
mass  directly  reduces  miterials  cost  and  in  general  improves  the  molr  H ty  of 
the  shielded  vehicle,  but  for  other  reasons  as  well.  For  instance,  ayaring 
also  improves  defense  against  projectiles  from  conventional  weapons  i  :  is 
well  known  that  a  single  sheet  of  steel  is  more  easily  perforated  tn  in  two 
separate  sheets  each  of  half  thickness)  and  lends  itself  well  to  tn  construc¬ 
tion  of  composite  shields  for  protection  against  mixed  gamma/neutroi  fields. 
Also,  as  will  be  shown  under  "Results",  the  shield  geometry  not  on).,  .effects 
the  total  dose  in  the  shielded  volume  but  also  the  dose  distribution  It  thus 
may  be  possible  to  create  local  "cold  spots"  at  regions  occupied  by  personnel 
merely  by  proper  1  nfiguration  of  the  shield,  not  by  Increasing  shie' j  rmss. 

The  subject  of  extension  of  the  Phase  I  work  will  be  treated  in  detail  in 
a  Phase  II  proposal;  we  wish  here  only  to  allude  to  the  richness  of  this  field 
of  investigation.  It  has  been  shown  in  this  work  that  the  relat've  sizes  and 
positions  of  uniform  shield  layers  of  constant  thickness  can  s 1 gr if  cant  1 y 
affect  dose  in  the  target;  however,  the  use  of  more  complex  geomett  es,  ruch 
as  conic-section  or  curved-surface  layers,  may  very  well  prove  to  live  . 
greater  effect.  For  instance,  consider  the  situation  portrayed  In  tig.  1 1 -la 
where  the  shield  layer  is  a  portion  of  a  cone  nappe.  A  photon  inte- acting  in 
the  conic  layer  will  generally  have  less  material  to  penetrate,  and  thus 
greater  chance  for  escape,  than  a  similar  photon  in  the  plane-parallel  layer 
of  Fig.  Il-lb.  Conic  or  other  types  of  shaped  layers  thus  appear  to  offer  an 
enhanced  potential  for  leakage  from  the  shield  after  scatter.  In  addition, 
density  and  composition  gradients  may  also  offer  advantages. 

For  the  protection  of  ground-vehicles,  we  note  that  a  series  of  thin 
curved  layers  as  shown  In  Fig.  1 1 -2a  might  afford  better  protection  tnar  one 
thick  shield  (Fig.  1 1 -2b ) .  A  photon  scattering  In  the  outer  layer  can  .tream 
through  the  annular  region  between  layers  and  either  escape  through  the  one 
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thin  layer  at  C  or  scatter  harmlessly  down  by  virtue  of  its  displacement  to 
point  C.  More  exotic  designs  are  possible  and  may  afford  increased 
protection. 


Finally  we  note  that  the  fact  that  the  MC  method  has  been  used  to  study 
the  problem  (and  in  fact  is  the  traditional  way  of  performing  detailed  dose 
calculations  for  armored  vehicles)  leads  to  the  interesting  possibility  of 
applying  the  Inverse  Monte  Carlo  (IMC)  method  [5]  to  the  design  optimization 
of  the  shields.  This  method,  developed  by  the  Principal  Investigator,  has 
been  applied  to  radiative  transfer  inverse  problems  [6]  and  to  the  design  of 
photon  filters  for  radiation  therapy  [7],  It  allows  complex  inve  >e  and 
optimization  problems  to  be  solved  In  a  noniterative  simulation  and  thus  has 
potential  application  to  efficient  gamma-ray  and  neutron  shield  design. 

B.  Geometry  Conventions 

The  geometry  we  consider  is  quite  simple.  We  assume  a  target  volume  Is 
suspended  in  free  space  and  subject  to  incident  gamma  radiation  which  origi¬ 
nates  at  a  distance  P  from  the  target.  For  finite  P  this  leads  to  a  slightly 
divergent  gamma-ray  beam;  as  P  becomes  infinite  this  reduces  to  a  uniform, 
parallel  beam.  The  assumption  of  a  vacuous  environment  reduces  the  number  of 
parameters  and  allows  the  feasibility  of  the  layering  principle  to  be  studied 
in  a  more  focused  and  isolated  sense. 

For  generality  and  completeness,  we  consider  the  target  to  consist  of 
either  a  right  circular  cylinder  or  a  sphere  and  score  the  dose  in  each 
simultaneously,  assuming  the  absence  of  the  other.  By  considering  a  rela¬ 
tively  large  cylinder  and  a  fairly  small  sphere,  we  can  construct  dose 
estimates  typical  of,  say,  a  cockpit  volume  and  a  central-axis  point.  Both 
targets  are  assumed  to  be  composed  of  water,  which  is  approximately  tissue- 
equivalent.  The  target  cylinder  has  radius  Rc  and  length  or  thickness  t;  its 
central  axis  is  col  inear  with  the  x-axis.  The  target  sphere  has  radius  Rs  and 
its  center  coincides  with  the  center  of  the  cylinder.  The  general  target- 
shield  geometry  is  shown  in  Fig.  1 1 - 3a . 

The  shield  is  composed  of  n  z  0  right  circular  cylindrical  layers  whose 
central  axes  are  collnear  with  the  x-axis  and  thus  with  the  central  axis  of 
the  target  cylinder.  The  layers  are  all  of  the  same  radius,  R^,  but  have 
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arbitrary  thicknesses  and  separations.  The  thickness  of  the  l^h  layer  is 
denoted  t^  and  the  separation  from  the  i-lst  member  is  denoted  s -j ,  as  shown  in 
Fig.  II-3a.  For  tMs  Phase  I  study,  the  shield  layers  consist  of  either  lead 
(high  Z,  high  densi:y)  or  water  {low  Z,  low  density)  and  in  a  given  case  all 
layers  must  be  of  tue  same  material  and  density. 

We  let  d  be  the  separation  between  the  target  cylinder  and  the  nearest 
face  of  the  closest  layer,  0  the  distance  from  the  target  cylinder  to  the 
front  face  of  the  fi*'st  layer,  T  the  total  shield  thickness,  S  the  total 
spacing,  X^,  the  x-coordinate  of  the  front  face  of  the  target  cylinder  and 
(Xs,0,0)  the  coordinates  of  the  center  of  the  target  sphere.  We  use  the 
convention  that  X^i  aid  X^2  specify  the  x-coordinates  of  the  front  and  back, 
respectively,  of  the  ith  shield  layer.  Then 


n 


T  =  2  ti  , 
i=l 

(i) 

n 

s  =  2  si  , 

i=l 

(2) 

D  =  Xt-si==Xt- 

xn  > 

(3) 

Xt  =  S  +  T  +  d  , 

(4) 

i  1 

xii  =  2  sj  -  2 
j=i 

j 

(5) 

and 

Xi2  =  Xn  +  t  .  (6) 

We  restrict  o.' selves  to  the  case  that  Rj,  =  Rc,  i.e.,  the  shield  layers 
have  the  same  radius  as  the  target  cylinder.  This  is  done  because  a  shield  of 
smaller  radius  woul  '  allov  direct  Irradiation  of  part  of  the  target  cylinder 
with  no  attenuation  and  a  shield  of  larger  radius  would  divert  some  photons 
toward  the  target  U  at  would  otherwise  stream  by  unaffected.  In  this  general 
geometry,  a  given  shield  configuration  is  described  by  the  number  of  layers, 
n,  the  target-to-shi  >ld  d 1  tance,  d,  and  the  individual  spacings,  Si,  and 
thicknesses,  t-j ,  for  a  tot.  1  of  2n  +  2  parameters. 


In  order  to  reduce  the  complexity  of  representing  a  shield  configuration, 
we  consider  the  simplified  geometry  of  Fig.  II-3b  which  is  subject  to  the 
restrictions  that  the  thicknesses  are  even, 

ti  =  T/n  ,  i  =  1,2,. ...n  ,  (7) 

and  the  spaclngs  are  even  except  for  the  first  which  Is  zero. 

si  =  $/(n  -1)  ,  1  =  2,3,. ..,n  (8a) 

=0  ,1=1  .  (8b) 

Under  these  restrictions,  we  can  represent  a  shield  configuration  by  the  four 
variables  T,  D  and  d.  This  geometry,  depicted  in  Fig.  1 1 -3b ,  is  referred 
to  as  the  simplified  geometry  and  will  be  used  to  represent  most  configura¬ 
tions  discussed  In  the  "Results"  chapter. 

It  is  noted  that  the  azimuthal  symmetry  of  the  geometry  considered  would 
allow  treatment  of  photon  transport  In  a  plane  only.  However,  extension  of 
this  work  to  more  complex  geometries  or  to  asymmetric  Incident  gamma-ray 
distributions  would  require  full  3-dimenslonal  (3-D)  transport  simulation. 
Thus,  the  analysis  and  code  development  has  used  a  3-D  transport  model. 
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III.  THE  MONTE  CARLO  MODEL 


The  direct  Monte  Carlo  method  provide  a  means  of  estimating  expected 
values  and  hence  definite  Integrals.  The  method  is  well  established  and  has 
been  used  extensively  in  radiation  transport  calculations  [e.g.,  8-12].  Thus, 
the  discussion  here  will  concentrate  on  its  application  to  the  problem  of 
interest  rather  than  on  the  basis  of  the  method  itself.  We  do  note,  however, 
a  few  general  points  in  order  to  establish  our  notation  and  conventions. 

The  Monte  Carlo  method  is  based  on  the  Strong  Law  of  Large  Numbers,  which 
basically  states  that  the  arithmetic  average  of  a  number  of  trials, 

.  1  N 

z  =  “  I  z(*i)  (1) 

N  1=1 

converges  almost  always  to  the  true  mean 
b 

<z>  =  /  z(x)  f(x)  dx  (2) 

a 

in  the  limit  of  large  N,  i.e., 

lim  z  =  <z>  ,  (3) 

N-to 

where  the  £i  are  sampled  from  the  probability  density  function  (pdf),  f(x). 
Furthermore,  the  Central  Limit  Theorem  provides  a  means  to  estimate  the 
variance  in  the  Monte  Carlo  estimator  z.  Much  of  the  power  of  the  Monte  Carlo 
method  derives  from  the  fact  that  variance  reduction  can  be  achieved  by 
various  means,  such  as  biased  sampling  and  contrived  estimators.  We  will 
describe  various  variance-reduction  procedures  as  they  are  applied  to  the 
layered  shield  problem,  in  particular  the  use  of  importance  sampling,  analytic 
equivalence  and  last-flight  estimators. 

The  geometry  we  will  consider  was  discussed  in  Chapter  II  and  is  shown  in 
Fig.  1 1 -3 .  In  the  following  sections  we  will  consider  photon  transport 
through  the  shield  layers  and  dose  deposition  in  the  targets. 
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A.  Photon  Transport 

First,  we  specify  the  initial  conditions.  These  consist  of  the  energy, 

E;  direction,  Q  =  (o,B,r)  with  a,  B,  y  the  direction  cosines;  position, 
(x,y,z);  and  statistical  history  weight,  W. 

Initial  Energy  and  Weight 

The  code  allows  two  options  for  specifying  Initial  energy:  mono- 
energetic  (ne  =  0)  and  histogram  (ne  >0).  If  ne  =  0,  the  source  energy  E0  Is 
input.  If  ne  >  0,  the  spectrum  cut-points,  Ej,  and  cumulative  fraction,  Sj, 
are  input  (see  Fig.  III-l).  The  energy  is  initialized  to 

E  *  E0  ,  ne  =  0  (4a) 

E  *  Ej_i  +  £(Ej  -  Ej-i )  ,  ne  #  0  ,  (4b) 

where  here  and  in  what  follows  £  is  a  random  variate  sampled  uniformly  from 
the  unit  Interval  (0,1).  The  particular  j  used  in  Eq.  4b  is  determined  by  the 
condition 

Sj— 1  <  £  <  Sj  ,  (5) 

where  S0  =  0  and  C  is  a  separate  uniform  random  variate.  The  history  weight 
is  initialized  to  unity, 

W  =  1  .  (6) 

Initial  Position  and  Direction 


The  Initial  position  is  determined  by  sampling  an  r  and  <(>  (see  Fig. 
Ill -2a )  so  as  to  give  a  uniform  distribution  over  the  area  of  a  circle,  using 

r  =  R,(5!)l/2  (7a) 


4>  -  2h£2  ’ 
and  then  letting 
x  =  Xn 
y  =  r  s1n$ 


z  =  r  cos<() 


Vi  .  A.1.  -  i  .  -  .  >  . .  j  -  „■  A.,  -  .  •  .  •  .  *  .  -  .  .  t .  «  . -  .  •  ■ 
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(a)  Front  View,  Showing  First  and  Last 
Annular  Regions  (of  Equal  Areas) 


Y 


Fi  rst 

Shield 

Layer 


(b)  Side  View,  Showing  Remote  Source  Point 


Figure  1 1 1 - 2 .  Geometry  of  Photon  Incident  on  First  Layer 


For  later  purposes,  we  wish  to  keep  track  of  the  radial  position  that  the 
photon  entered  the  first  layer.  To  accomplish  this,  we  construct  annuli  of 
equal  area  on  the  front  face  of  the  first  layer,  as  in  Fig.  Ill -2a ,  and  store 
the  annulus  number,  I,  corresponding  to  the  entry  point  such  that 
rj_i  <  r  <  rj,  where  r i  is  the  outer  radius  of  annulus  i.  These  radii  are 
determined  by  requiring  that  the  area  of  each  annulus  be 


Ao  = 


"  R*2 


(9) 


(10) 


where  nd  is  the  number  of  annuli.  Then, 

H  =  [rj2  +  n-!?]1/2  ,  i  >  1  , 

with  r0  =  0. 

The  initial  direction  is  determined  by  assuming  the  radiation  originates 
at  a  point  a  stance  P  from  the  first  shield  layer  (Fig.  Ill -2b ) .  Under 
these  conditions,  the  initial  direction  cosines  are  given  by 

a  =  cose  (11a) 

&  =  sine  si nm  (lib) 

y  =  sine  cos<t>  (11c) 

where 

9  =  tan"l  (r/P)  ,  (12) 

and  r  and  tp  are  given  by  Eqs.  7.  In  the  limit  of  infinite  P,  the  beam  is 
parallel  and  Eqs.  11  reduce  to  a  =  1 ,  8  =  y  =  0;  this  is  the  normal  case  we 
have  treated  although  Eqs.  11  were  used  in  one  study  to  consider  a  diverging 
beam. 


Oetermi ne  Position  and  Type  of  Interaction 

In  order  to  determine  how  far  the  photon  travels  before  interacting, 
we  must  determine  the  total  linear  attenuation  coefficient,  of  the  shield 
material.  (In  this  Phase  I  study  we  have  assumed  that  all  shield  layers 
consist  of  the  same  materials;  however,  the  code  has  been  written  so  that 
different  compositions  and  densities  can  be  input  for  the  various  layers  so 
that  in  Phase  II  more  elaborate  shield  conf iqurations,  for  instance  having 


atomic  number  or  density  gradients,  can  be  studied.)  We  will  also  need  the 
various  partial  cross  sections  for  Compton  (incoherent)  scattering,  a,  photo¬ 
electric  absorption,  t,  and  pair  production,  in  order  to  determine  the  type 
of  interaction.  Thus,  we  call  a  subroutine  that  calculates,  using  piecewise 
polynomial  approximation,  the  cross  sections  a,  x,  *•  and  their  sum 

p(E)  =  a  +  x  +  r  ,  (13) 

given  the  photon  energy,  £.  We  have  not  included  Rayleigh  (coherent)  scat¬ 
tering  nor  bound  electron  effects  on  the  Compton  process  in  this  Phase  I  study 
because  these  involve  small-angle  scattering  and  the  Rayleigh  cross  section  is 
generally  nail;  hence  these  effects  are  not  needed  in  order  to  establish 
feasibllu^.  However,  both  effects  would  be  included  in  the  Phase  II  study 
since  optimizing  shield  design  using  the  layering  principle  would  require 
proper  treatment  of  even  small-angle  scatters. 

The  distance  to  interaction,  s,  is  obtained  by  sampling  from  the  pdf 


f(s)  =  v  e-Ps  ,  (1A) 

from  which  we  obtain 
Jin  i 

s  - - •  (15) 

v- 

The  coordinates  of  the  photon  are  easily  updated  to  the  new  Interaction  point 
using 

x'  =  x  +  $a  (16a) 

y'  =  y  +  sb  (16b) 

z'  =  z  +  sr  •  (16c) 


In  applying  Eqs.  16,  we  check  to  see  that  s  is  less  than  the  distance,  A,  to 
the  nearest  boundary  (see  Fig.  1 1 1 - 3 ) .  If  so,  Eqs.  16  are  used  directly;  if 
not,  the  photon  has  escaped  the  shield  layer  and  we  must  perform  t  e  following 
procedure: 

(1)  If  the  photon  exits  through  the  cylindrical  surface,  the 
history  is  terminated  (path  1  in  Fig.  1 1 1 - 3 ) . 


(2)  If  the  photon  escapes  through  a  planar  surface,  we  update  the 
photon  position  to  the  exit  point,  reduce  the  distance  s  by  the 
distance  through  this  layer  A,  find  the  intersection  with  the 
next  layer  (if  any),  position  the  photon  at  this  entrance 
point,  apply  Eqs.  16  with  (x,y,z)  representing  the  entrance 
point  and  s  the  path  length  reduced  by  the  distance  A,  and 
repeat  the  process  (path  2  in  Fig.  III-3).  If  there  is  no 
intersection  with  another  layer,  the  history  is  terminated 
(path  3  in  Fig.  1 1 1-3 ) . 

The  photon  paths  from  incidence  on  the  first  shield  layer  to  first 
interaction,  from  one  interaction  to  another  or  from  one  Interaction  to  entry 
into  one  of  the  targets  are  straight-line  segments.  To  determine  where  a 
given  pat!  ntersects  shield  or  target  boundaries,  we  must  be  able  to  find  the 
intersection  of  a  ray  starting  at  (x,y, z)  and  having  direction  (a,B,y)  with 
planar,  cylindrical  and  spherical  surfaces.  We  let  the  Intersection  point  be 
(xc,yc,zc)  and  occur  at  radial  distance  rc  from  the  x-axis.  The  determination 
of  these  intersection  points  is  straightforward;  for  completeness  the  solu¬ 
tions  are  presented  in  Appendix  A.  The  distance  between  (x,y,z)  and  (xc,yc, 
zc)  is  then  simply 

4  *  t(xc  -  *)2  +  (yc  -  y)2  +  <*c  -  z)2]172  •  (in 

Force  Scatter 


Since  the  pair  production  interaction  is  inevitably  followed  by 
annihilation  of  the  positron  and  emission  of  two  0.511  MeV  gamma-rays,  we 
treat  it  for  shielding  purposes  as  a  scattering  interaction.  We  refer  to  this 
"scatter"  as  pair  production/annihilation  (PP/A).  This  Interaction  leads  to 
one  0.511  MeV  photon  emitted  Isotropically  and  a  second  0.511  MeV  photon 
antiparallel  to  the  first.  The  other  scattering  interaction  we  treat  is  the 
incoherent  scatter,  referred  to  as  Compton  (C).  We  do  not  include  the 
scatter-function  correction  for  bound  electron  effects  and  we  neglect  coherent 
(Rayleigh)  scattering.  The  final  Interaction  we  consider  Is  the  photoelectric 
effect  (P),  which  is  the  complete  transfer  of  the  photon  energy  to  a  bound 
electron. 


We  force  scatters  as  a  variance  reduction  measure  to  Improve  the  effi¬ 
ciency  of  the  simulation.  Since  P  interactions  remove  photons  from  considera¬ 
tion,  they  can  no  longer  deposit  energy  in  the  target,  the  event  we  are 
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as 


interested  in  scoring.  Thus,  we  pick  PP/A  interactions  with  frequency 
r/(x  +  tr)  and  C  scatters  the  rest  of  the  time.  This  requires  that  we  weight 
the  history  by  the  factor 

t  +  tr 


which  is  the  probability  of  scatter.  All  cross  sections  are  of  course  deter¬ 
mined  at  the  current  photon  energy. 

Find  Post-Scatter  Conditions 

If  £  <  r/(t  +  r),  we  treat  the  interaction  as  an  Isotropic  scatter. 
Letting  <|>  be  azimuthal  angle  and 

a  =  cos  e  (19) 


where  e  is  the  scatter  angle,  we  sample  to  obtain 
-  =  1  -  2  d  , 


(20a) 

(20b) 


<t>  =  2*  $2  .  (20b) 

then  update  the  direction  cosines  using  the  formulas  of  Selph  and  Garrett  [11] 


sin  e 

a'  =  a  to  +  -  [a  y  C0S<)>  -  B  sin$] 

r 

sin  e 

B '  =  B  to  +  -  [B  y  cos<J>  +  a  s i n«J>  j 

r 

y'  =  y  (u  -  T  sine  cos<(»  , 


(21a) 


(21b) 


(21c) 


where 


T  =  (1  -  r2)1/2  (2i 

and  (ot',B',y')  are  the  direction  cosines  after  the  scatter.  We  save  the 
negative  values  of  these  direction  cosines  and  the  values  of  W,  x,  y,  and  z 
for  later  simulation  of  the  second  photon  resulting  from  the  positron 
annihi lation. 


•■'.wNKs 


If  £  >  ic/(t  +  *-),  we  treat  the  interaction  as  a  Compton  scatter.  We 
sample  the  azimuthal  angle  as  In  Eq.  20b  and  determine  the  cosine  of  the 
scatter  angle  in  the  following  manner.  We  sample  «  from  the  normalized 
version  of  g(«|n)  in  Eq.  8-9b  (see  Appendix  B), 

9(«ln) 

9n(«ln)  =  -  ,  (23) 

0 

where  D  is  given  by  Eq.  B-14,  and  use  Eq.  B-5  to  obtain  to,  i.e., 
n  +  1  -  a 

«  =  -  ,  (24) 

n 

where 

E 

n  =  -  .  (25) 

0.511 


The  procedure  used  to  sample  from  gn  is  a  version  of  the  Kahn  [8]  rejection 
technique  taken  from  Dunn  and  Gardner  [13].  We  then  update  the  direction 
cosines  according  to  Eqs.  21  and  determine  the  post-scatter  photon  energy  from 

E*  =  E/«  .  (26) 

Whether  a  PP/N  or  C  scatter  was  chosen,  we  test  to  see  if  the  product  of 
history  weight  and  photon  energy  is  larger  than  some  input  cutoff  value,  e. 

If  not,  we  terminate  the  photon  history  since  the  additional  contribution  to 
the  dose  in  the  target  would  be  negligible;  if  so,  we  continue  the  simulation 
by  calculating  the  cross  sections  at  the  new  energy  and  proceeding  as  dis¬ 
cussed  in  the  section  "Determine  Position  and  Type  of  Interaction." 

B.  Target  Scoring 

We  score  two  target  dose  quantities,  one  typical  of  the  central -axis, 
first-collision  dose  and  one  representatl ve  of  the  total  dose  deposited  in  a 
larger  mass.  The  descriptors  are  actually  fractional  energy  depositions, 
representing  energy  deposited  In  the  target  per  unit  energy  Incident  on  the 
shield  front  face.  The  central-axis  descriptor  is  defined  as  the  expected 
energy  deposited  in  the  target  sphere  by  first  collisions  in  the  sphere  per 
unit  energy  incident  on  the  front  shield  face;  it  is  denoted  Fs  and  is 
expressed  In  units  of  eV/MeV  (eV  deposited  per  MeV  Incident).  By  first- 
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collision  we  mean  that  any  secondary  photons  born  in  the  first  collision 
within  the  sphere  are  assumed  to  escape  the  sphere.  On  the  other  hand,  all 
energy  given  to  electrons  is  assumed  to  remain  in  the  sphere.  The  quantity 
scored  for  the  larger  cylindrical  target  is  denoted  Fc  and  is  defined  as  the 
total  energy  deposited  in  the  cylinder  (in  keV),  including  that  from  secondary 
photons,  per  MeV  incident  on  the  first  shield  layer.  Thus,  for  the  target 
cylinder  we  track  both  primary  and  secondary  photon  transport  through  the 
cylinder. 

We  score  the  quantity  Fc  using  a  last-flight  estimator.  Thus,  for  every 
photon  path,  including  the  incident  one,  we  extend  the  path  to  see  if  it 
intersects  the  target  cylinder.  If  so,  we  determine  the  total  distance,  «, 
through  all  shield  layers  along  this  path  and  allow  the  photon  to  enter  the 
cylinder  with  weight  e“P5,  where  is  the  total  linear  attenuation  coefficient 
of  the  shield  material  at  the  current  photon  energy.  We  then  force  an  inter¬ 
action  in  the  target  with  weight  1  -  e-^6^  where  pt  is  the  total  linear 
attenuation  coefficient  of  the  target  material  at  the  current  photon  energy 
and  «£  is  the  distance  through  the  target  to  escape. 

At  each  interaction  point  in  the  target  cylinder  we  score  average  energy 
deposited  with  the  estimator  [a^T  +  hE  +  *^(£-1 .022) ]/nt>  where  T  Is  the 
average  energy  deposited  in  a  Compton  collision  at  the  current  photon  energy, 
E,  and  the  subscript  t  refers  to  the  target  cylinder.  The  determination  of  T 
is  discussed  in  Appendix  B.  We  then  force  a  PP/A  or  C  interaction,  with 
weight  (»t  +  sample  to  determine  which  of  the  two  occurs  and  continue 

simulation  through  the  target,  scoring  average  energy  deposition  at  each 
interaction  point  and  terminating  when  the  product  of  photon  weight  and  energy 
is  reduced  below  the  cutoff,  e.  After  termination,  we  return  to  the  original 
interaction  point  in  the  shield  and  resume  photon  transport  through  the 
shield.  Thus,  the  cylinder  energy  deposition  estimator  for  a  given  history 
has  the  form 


-u«  v  ((n  +  -ut«t  fatT  +  TtE  +  ^(E-1-022)]) 

AEC  =  e  **  l  -  (1  -  e  “t«t)  - L  ,  (27) 

j=i  Mt.  Mt  r 


3~y  {  Mt  Pt 

where  v  is  the  number  of  interactions  in  the  target  before  history  termination 


We  score  the  central-axis  quantity  Fs  in  a  slightly  different  manner. 
Because  the  sphere  is  small  the  probability  that  a  photon  path  will  intersect 
the  sphere  is  small.  To  improve  the  simulation  efficiency,  we  find  the  solid 
angle  subtended  by  the  sphere  at  each  interaction  point  in  the  shield  and 
force  scatters  into  this  solid  angle.  This  procedure  is  discussed  in  Appendix 
C.  We  then  extend  the  path  to  the  sphere,  force  an  interaction  there,  and 
score  the  energy  deposition  with  an  estimator  like  that  used  for  the  cylinder. 
However,  we  do  not  continue  the  photon  simulation  through  multiple  interac¬ 
tions  in  the  sphere,  because  of  its  small  size.  Rather,  we  return  to  the 
interaction  point  in  the  shield  and  continue  simulation  there.  The  sphere 
energy  deposition  estimator  thus  has  the  form 


where  We  and  W^  are  weight  factors  for  forcing  scatter  toward  the  sphere  and 
are  defined  in  Appendix  C. 


We  have  not  attempted  to  develop  highly  sophisticated  dose  descriptors 
that  require  use  of  materials  more  tissue-equivalent  than  water  (for  instance, 
as  specified  by  the  ICRU  [14])  or  that  are  truly  dose-equivalent  [15]  or 
representati ve  of  actual  dosimeter  measurements,  since  these  are  unnecessary 
to  study  the  gross  effects  of  shield  layering  and  placement.  However,  in 
studying  the  finer  details  of  shield  design,  a  more  sophisticated  target 
scoring  scheme  could  be  employed. 
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IV.  CODE  VERIFICATION 


A.  General 

We  have  verified  the  components  of  the  code  and  Its  internal  consistency 
in  various  ways.  In  general,  the  code  was  written  in  modular  form  and  each 
module  separately  tested  and  checked  carefully  before  integration  with  the 
rest  of  the  code.  For  instance,  individual  subroutines  were  written  to  find 
the  intersections  of  rays  with  planes,  cylindrical  surfaces  and  spheres  (as 
described  In  Appendix  A)  and  each  subroutine  was  tested  before  inclu'  ''n  In 
the  code.  The  subroutine  to  determine  conditions  following  Incoherent  scatter 
was  verified  by  calling  the  subroutine  repeatedly  for  the  same  initial  photon 
energy,  plotting  the  frequency  distribution  of  the  resulting  scattering  angle 
and  comparing  to  the  known  distribution  given  by  Evans  [16].  The  photon  cross 
section  library  was  constructed  by  piecewise  polynomial  fit  over  the  range  10 
keV  to  20  MeV  to  the  data  of  Hubbell  [17].  It  was  thoroughly  tested  and  the 
relative  error  between  model  and  data  was  shown  to  be  generally  less  than  3 
percent. 

The  Integrated  code  was  also  tested.  We  have  shown  that  gamma-ray 
transport  through  two  contiguous  shield  layers  is  Identical  to  transport 
through  one  layer  of  the  same  total  thickness,  as  a  partial  check  of  the 
geometry  of  treating  layers.  We  also  ran  several  MC  histories  with  a  detailed 
output  of  various  intermediate  variables  which  were  checked  Independently 
using  a  hand-calculator. 

As  a  further  check  of  the  code,  we  have  generated  the  various  surface 
quantities  required  In  the  Inverse  formulation  of  Siewert  and  Dunn  [18]  for  a 
monoenergetic,  monodirectlonal  beam  Incident  at  a  point  on  the  surface  of  a 
single  shield  layer  of  radius  1000  cm  and  thickness  4  cm,  approximating  a 
plane  parallel  medium  of  Infinite  radius.  Since  this  Inverse  solution  assumes 
Isotropic  scattering,  we  used  the  isotropic  scattering  option  available  in  the 
code,  which  keeps  the  Compton  cross  section  at  its  value  at  the  incident 
gairma-ray  energy  and  samples  the  cosine  of  the  scatter  angle  from  a  uniform 
distribution  on  (-1,1).  For  this  case,  one  form  of  the  Inverse  solution  can 
be  expressed  [18] 


29 
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(1) 


a  8u  ii>(-a0,<j>0) 
p 

where  a  and  p  are  the  Compton  and  total  linear  attenuation  coefficients, 
respectively,  o0  >  0  is  the  x-direction  cosine  of  the  Incident  beam,  $0  is  the 
azimuthal  angle  of  the  incident  beam  measured  positive  CCW  from  the  +z  axis, 
^(-"o^o)  is  ex1t  angular  flux  in  the  direction  opposite  the  incident  beam 
integrated  over  the  entire  left  surface,  ity  is  the  total  flux  integrated  over 
all  directions  over  the  entire  left  surface,  and  i|»r  is  the  total  flux  Inte¬ 
grated  ove>  all  directions  over  the  right  surface.  The  required  exit  fluxes 
are  found  by  a  last  flight  estimator,  where  the  estimator  after  each  inter¬ 
action  (including  initial  entry  to  the  shield)  has  the  form 

W  e~PA 

♦  =  -  >  (2) 

id 

where  W  is  the  incoming  history  weight,  A  the  distance  to  the  shield  surface 
and  c  the  cosine  of  the  angle  between  the  extended  gamma-ray  path  and  the 
outward  normal  to  the  shield  surface  at  the  exit  point.  Substituting  the 
simulated  quantities  into  Eq.  1,  we  have  recovered  the  ratio  of  scattering  to 
total  cross  sections  at  the  energy  of  the  incident  gamma-rays  to  within  one 
percent  of  the  correct  value.  This  check  verifies  the  geometrical  treatment 
of  photon  transport  and  the  implementation  of  last-flight  estimation. 
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V.  RESULTS 


For  our  analysis  we  have  chosen  a  standard  target  arrangement  consisting 
of  a  one-meter  diameter  by  one-meter  long  circular  cylinder  and  a  5-cm 
diameter  sphere  located  at  the  center  of  the  cylinder.  Both  targets  are 
assumed  to  consist  of  water,  which  is  approximately  tissue-equivalent.  The 
cylinder  is  representative  of  a  relatively  large  target  within  which  multiple 
photon  interactions  may  occur;  the  sphere  is  small  enough  that  multiple  photon 
interactions  inside  it  are  unlikely  and  that  its  dose  is  a  reasonable  measure 
of  central-axis  dose.  The  fractional  energy  absorbed  per  incident  MeV,  as 
discussed  in  Chapter  III,  is  determined  for  each  target,  in  the  absence  of  the 
other. 

The  assumed  photon  spectrum  is  taken  from  Glasstone  and  Dolan  [19], 
pp.  360-1,  and  is  approximated  by:  Eq  =  0,  E}  =  0.75,  E2  =  2,  E3  =  4.5, 

E4  =8,  E5  =  12;  Sj  =  0.7,  $2  =  0.8,  S3  =  0.89,  S4  =  0.975,  S5  =  1.  We  have 
chosen  to  treat  three  incident  energy  conditions:  monoenergetic  low,  mono- 
energetic  high  and  the  full  Glasstone  and  Dolan  [19]  spectrum.  In  this  way 
low-  and  high-energy  effects  can  be  isolated  initially  for  separate  study  and 
then  considered  in  a  composite  manner  at  a  later  stage.  The  midpoint  of  the 
first  energy  interval,  which  contains  70  percent  of  the  photons,  is  0.375  MeV 
and  so  defines  the  low-energy  case.  The  monoenergetic  high  case  is  taken  as 
2.5  MeV,  which  is  the  approximate  weighted  average  of  the  remaining  30  percent 
of  the  spectrum. 

It  is  well  known  that  gamma  radiation  in  a  thin-beam  geometry  is  atten¬ 
uated  exponential ly.  In  broad-beam,  thick-shield  geometry,  however,  bu41d-up 
occurs  due  to  multiple  scattering  and  the  dose  (or  intensity)  on  the  far  side 
of  the  shield  is  larger  than  a  strict  exponential  model  of  the  form 

0(T)  =  0o  e'^T  (1) 

would  predict.  In  Eq.  1  D(T)  denotes  some  dose  measure  at  the  far  surface  of 
the  shield  of  thickness  T,  D0  the  same  measure  for  T  =  0  and  y  the  linear 
attenuation  coefficient  of  the  shield  material.  The  ratio  of  actual  dose  to 
dose  estimated  by  the  simplified  model  of  Eq.  1  Is  commonly  called  the  dose 
buildup  factor.  This  factor  depends  on  shield  material,  photon  energy  and 


shield  thickness,  but  can  b*.  ^s  large  as  10  or  20  or  even  100  [4,20]  and  tends 
to  increase  with  increasing  shield  thickness.  The  fact  that  this  buildup  of 
dose  can  be  so  significant  and  that  it  is  due  to  scattering  in  the  shield  beg 
investigation  of  the  extent  to  which  the  scattered  dose  contribution  can  be 
reduced  by  modifying  either  the  scattering  process  or  the  transport  of  photons 
between  interactions.  Various  aspects  of  shield  design  are  studied  in  the 
following  sections. 

Effect  of  Shield  Thickness 

In  Table  V-l  we  show  the  quantities  Fs  and  Fc  as  functions  of  shield 
thickness  for  a  shield  composed  of  a  single  layer  of  lead  positioned  at 
various  target-to-shield  distances,  d.  The  central  sphere  quantities,  Fs,  are 
plotted  in  Fig.  V-l.  Several  things  are  apparent  from  these  results.  First 
of  all,  the  variation  of  Fs  with  thickness  is  approximately  exponential  at 
large  thicknesses  (straight  lines  on  the  semi logari thmic  scale),  especially 
for  the  low-energy  cases  (Cases  1,  2  and  3).  However,  for  small  thicknesses, 
the  behavior  is  decidely  not  exponentially  decreasing,  especially  for  d  =  0 
(i.e.,  for  the  shield  near  the  target),  and  even  for  large  thicknesses  the 
rate  of  decrease  Is  much  smaller  than  the  attenuation  rate  of  primary  photons, 
e-^,  where  p  =  2.84  cm"!  at  0.375  MeV  and  n  =  0.044  at  2.5  MeV  (Cases  a  and 
b  in  Fig.  V-l).  In  fact,  at  2.5  MeV  there  is  actually  an  Increase  in  target 
dose  for  thin  shields  over  the  no-shield  case,  clearly  pointing  out  the 
scattered-photon  contribution. 

From  Table  V-l  it  Is  clear  that  the  cylinder  dose  quantity,  Fc,  behaves 
similarly,  although  there  is  no  increase  in  dose  over  the  no-shield  case  for 
thin  shields,  even  at  2.5  MeV.  It  is  also  obvious  from  these  results  that  the 
farther  the  shield  Is  from  the  target,  the  lower  the  dose.  This  effect  will 
be  investigated  next. 

Effect  of  Target-to-Shield  Distance 

We  now  consider  a  single  shield  of  fixed  thickness,  T,  and  study  the 
effect  on  dose  of  where  the  shield  is  placed  relative  to  the  target.  Several 
cases  are  considered  in  Table  V-2,and  Fig.  V- 2.  At  T  =  2  mfp  and 
E0  =  0.375  MeV  (Case  5)  we  see  that  a  reduction  in  the  quantity  Fs  by  almost  a 
factor  of  three  is  possible  simply  by  separating  the  shield  from  the  cylinder 


TABLE  V-l.  EFFECT  OF  SHIELD  THICKNESS  AT  VARIOUS 
TARGET-TO-SHIELD  DISTANCES,  d,  FOR  LEAD-SHIELD  LAYERS 


Figure  V-2.  Variation  of  Sphere  Energy 
Fraction  with  Target-to- Shi»M  Di  stance 


target  from  0  to  4  meters.  At  4  mfp  {Case  6)  the  reduction  is  greater  than  a 
factor  of  three.  This  effect  is  exaggerated  at  higher  energies  (in  Case  7 
the  reduction  is  a  factor  of  four),  presumably  because  of  the  approximate  iso¬ 
tropy  of  annihilation  interactions  as  opposed  to  the  forward-biased  scattering 
of  Compton  interactions. 

Case  8  In  Table  V-2  shows  that  while  the  absolute  doses  are  considerably 
larger  for  a  low-Z,  water  shield  of  2  mfp  thickness  the  effect  of  target-to- 
source  distance  is  significantly  increased,  with  the  ratio  Fs(d  =  0)/ 

Fs(d  =  4m)  being  greater  than  6.  This  is  presumably  due  to  the  fact  that 
hydrogen  is  a  more  effective  Compton  scatterer  than  other  elements,  since 
a  a  Z/A  and  Z/A  =0.5  for  all  elements  except  hydrogen  for  which  Z/A  =  1  (here 
Z  is  atomic  number  and  A  atomic  weight).  This  suggests  that  a  shield  that 
includes  hydrogen  (as  in  a  polymer)  at  large  D  to  induce  scattering  and  then  a 
high  density  material  (such  as  steel  or  lead)  to  provide  high  attenuation 
might  be  more  effective  than  a  strictly  high-Z  shield  of  equal  mass.  We  note 
also  that  such  a  conf iguration  Is  ideal  for  shielding  against  neutrons  (the 
hydrogen-rich  layer  would  thermal ize  the  neutrons  which  could  then  be  absorbed 
by  a  thin  layer  of  cadmium  or  of  borated  or  rare-earth  doped  material). 

Effect  of  Layering 

The  effect  of  layering  cannot  be  directly  Isolated  from  the  effect 
of  target-to-shield  separation,  since  any  splitting  of  a  shield  into  noncon¬ 
tiguous  layers  necessarily  changes  the  target-to-shield  separation  for  at 
least  one  of  the  layers.  Thus,  we  consider  various  aspects  of  shield  layering 
separately.  First,  we  show  In  Table  V-3  the  dependence  of  Fs  and  Fc  on 
spacing,  S,  for  a  shield  split  into  two  layers  with  the  first  layer  at 
constant  d.  Cases  9  and  10  (E0  =  0.375  MeV)  demonstrate  the  general  dose 
decrease  with  Increasing  S;  whether  the  first  layer  is  close  to  or  removed 
from  the  target,  there  is  an  approximately  20  percent  decrease  in  Fs  in  going 
from  $  =  0  (effectively,  one  thick  shield)  to  S  =  3  m.  These  resu-ts  indicate 
that  the  effectiveness  of  layering  tends  to  decrease  as  d  increases  but  that 
the  overall  shield  effectiveness  Improves.  Case  11  shows  that  the  effect  is 
greater  at  2.5  MeV  and  Cases  12  and  13  show  that  the  effect  Is  much  greater 
for  water  than  lead.  The  results  are  plotted  in  Fig.  V-3. 


37 


TABLE  V-3.  EFFECT  OF  SPLITTING  SHIELD  INTO  2  LAYERS  (n  =  2) 


Case 

Eo 

S 

T 

D 

d 

Es 

No. 

(MeV) 

(cm) 

(mfp) 

(cm) 

(cm) 

(eV/MeV) 

Lead  Shi  el di n 


0.375 


0.375 


0.375 


0 

98.59 

198.59 

298.59 

0 

98.59 

198.59 

298.59 

0 

191.74 

291.74 

0 

199.295 

399.295 


1.41 

0 

■n 

100 

0 

200 

0 

300 

0 

KH 

201.41 

300 

7.6 

6.8 

400 

200 

6.4 

500 

200 

6.2 

208.26 

200 

22.4 

400 

200 

15.9 

500 

200 

14.4 

400.705 

400 

35.2 

600 

400 

33.9 

800 

400 

33.4 

Water  Shieldln 


13  0.375  0 


2 

418.4 

400 

77.7 

2 

600 

400 

55.3 

2 

800 

400 

51.7 

The  effect  of  splitting  into  many  layers  is  Indicated  in  Table  V-4  and 

Fig.  V-4.  Case  14  Is  for  constant  spacing  between  layers,  S/(n  -  1),  while 

Cases  14  and  16  allow  the  spacing  between  layers  to  decrease  with  increasing 
n.  It  is  apparent  that  the  effect  is  again  greater  for  high  energy,  but 
significant  for  both  low  and  high.  We  note  that  the  reduction  with  n  would  be 
greater  in  Cases  15  and  16  if  S/(n  -  1)  were  constant  at  200  cm.  As  a 
contrived  means  to  study  the  effect  in  the  limit  as  n  -»  ®,  we  ran  a  case  in 

which  we  spread  an  hypothetical  lead  shield  over  a  2-m  distance  and  lowered 

the  density,  to  keep  the  mass  constant.  This  result  supports  the  conclusion 
that  for  fixed  d  and  0  (i.e.,  under  constraints  on  minimum  and  maximum 
target-to-shield  separations),  there  is  little  gain  in  splitting  into  more 
than  about  4  layers,  i.e.,  although  layering  always  helps,  its  marginal 
Improvement  is  small  after  about  n  =  4  layers. 

Effect  of  Thickness  Gradient 

We  have  shown  that  layering  (including  the  effects  of  varying  D)  can 
lead  to  substantial  dose  reduction  (e.g.,  greater  than  a  factor  of  2  in  Case 
15).  We  now  wish  to  Investigate  the  prospects  for  dose  reduction  under  the 
constraint,  likely  to  be  encountered  in  actual  shielding  situations,  that  both 
d  and  D  are  constrained,  i.e.,  the  distance  from  the  target  to  the  first 
(farthest)  shield  cannot  exceed  some  limiting  value.  The  results  of  this 
study  are  shown  in  Table  V-5,  and  are  quite  impressive.  Case  17  demonstrates 
that  at  0.375  MeV  an  11  percent  decrease  in  Fs  results  when  the  relative 
proportion  of  shield  thicknesses  Is  changed  from  1:1:2  to  2:1:1,  with  the 
shield  members  at  the  same  locations  In  the  two  cases.  At  2.5  MeV  (Case  18), 
this  decrease  Is  almost  25  percent  for  the  same  proportionate  thickness  change 
and  is  greater  than  42  percent  in  going  from  thickness  proportions  1:1:2  to 
6:2:1.  The  effect  Is  similarly  impressive  for  water  shields  (Case  19),  as 
expected.  Hence,  It  appears  that  rather  substantial  reductions  in  central- 
axis  dose  are  possible,  even  when  structural  or  other  constraints  dictate 
that  part  of  the  shield  must  be  positioned  close  to  the  target. 


V**  S 


T 

(mfp) 


201.41 

400 

800 

1200 

1600 

208.26 

400 

600 

1000 

1.41 

200 

200 

200 


A  hypothetical  lead  shield  of  thickness  200  cm,  density  0.079947  gm/cm3,  w 
modeled,  having  the  same  total  mass  (126.88  kg)  as  the  1.41  cm  shield  with 
normal  density  11.34  gm/cm3.  This  simulates  the  stretching  out  of  the 
shield  into  an  infinite  number  of  layers. 


d 

(cm) 

Fs 

(eV/MeV) 

■a 

7.6 

6.4 

200 

5.4 

200 

5.0 

200 

4.9 

200 

22.4 

200 

15.9 

200 

11.4 

200 

8.7 

0 

18.9 

0 

15.0 

0 

14.2 

0 

13.7 
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TABLE  V-5.  EFFECT  OF  THICKNESS  GRADIENT; 
d  =  200  cm,  D  *  800  cm,  n  =  3 


(keV/MeV) 


a.  Lead  Shielding,  T  =  1.41  (4  mf 


0.375 

.3525 

.3525 

.705 

6.1 

.282 

.564 

.564 

5.7 

.470 

.470 

.470 

5.6 

.705 

.3525 

.3525 

5.4 

2.5 

2.065 

2.065 

4.13 

13.3 

1.652 

3.304 

3.304 

11.9 

2.7534 

2.7533 

2.7533 

10.9 

4.13 

2.065 

2.065 

10.0 

5.509 

1.834 

0.917 

7.6 

Water  Shieldln 


=  18.4  cm 


1  19 

0.375 

4.6 

4.6 

9.2 

70.1 

6.134 

6.133 

6.133 

56.6 

9.2 

4.6 

4.6 

53.0 

Effects  for  the  Who le-Spectrum  Case 

The  results  presented  so  far  were  for  monoenergetic  incident 
photons.  We  view  typical  cases  in  both  the  simplified  and  generalized 
geometries  in  Table  V-6  for  the  Glasstone  and  Dolan  [19]  spectrum.  Here,  the 
sharp  decrease  of  Fs  and  Fc  with  increasing  source-to-closest-layer  distance, 
d,  is  evident  in  Case  20,  plotted  in  Fig.  V-5a,  and  the  substantial  decrease 
with  increasing  number  of  layers  (for  constant  d)  is  clear  from  Case  21,  shown 
in  Fig.  V-5b.  The  effect  of  reducing  the  thickness  of  the  layer  nearest  the 
target  to  reduce  dose  is  obvious  from  Case  22. 

Effects  of  a  Diverging  Incident  Beam 


The  results  of  Table  V-7  indicate  that  the  dose  reduction  effects 
are  similar  for  a  diverging  as  opposed  to  a  parallel  beam.  In  general,  the 
trends  are  similar  in  the  two  cases,  but  the  dose  quantities  for  the  diverging 
case  are  approximately  80-90  percent  of  those  for  the  parallel  case,  as 
expected.  The  specific  case  considered  involved  significant  divergence,  since 
the  source  was  assumed  to  be  only  100  m  from  the  front  of  the  cylindrical 
target.  Thus,  we  conclude  that  for  beams  typical  of  remote  bursts  or  directed 
energy  weapons  (which  should  be  bounded  by  the  parallel-  and  divergent-beam 
cases  considered  here),  the  effects  of  layering,  thickness  gradients  and 
shield  geometry  are  quite  significant. 


BLE  V-6.  RESULTS  FOR  INCIDENT  PHOTONS  DISTRIBUTED 
CCORDING  TO  THE  GLASSTONE  AND  OOLAN  1191  SPECTRUM 


(eV/M eV) 


122.4 

100.9 

66.9 

51.7 

44.2 


Effect  of  Thickness  Gradient,  d 


Geometry 


=  200  cm,  D  =  800  cm.  Generalized 


tl 

(cm) 

t2 

(cm) 

t3 

(cm) 

Fs 

(eV/MeV) 

1.6 

1.2 

1.2 

45.4 

1  1 

1 

1 

43.0 

2 

.5 

.5 

38.2 
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VI.  CONCLUSIONS 


This  Phase  I  study  draws  attention  to  the  importance  of  geometric 
optimization  in  shielding  both  terrestrial  and  space  vehicles  against  directed 
gamma-ray  beams.  By  means  of  a  versatile  and  efficient  Monte  Carlo  model 
developed  for  this  purpose,  it  has  been  shown  that  doses  to  the  target  may  be 
reduced  by  factors  of  2  or  3  without  increase  in  shield  mass,  through  careful 
choice  of  shield  configuration;  al ternati vely,  for  a  given  acceptable  dose, 
the  shield  mass  may  be  considerably  reduced.  The  tremendous  influence  of 
target-shield  geometry  has  been  demonstrated  and  the  effects  of  target-to- 
nearest-1.-  '  distance,  layer  spacing,  number  of  layers  and  layer  thickness 
gradient  nave  been  quantified.  The  studies  to  date  also  indicate  promising 
areas  for  further  improvements.  Composition  gradients,  wherein  light-weight, 
effective-scattering  materials  would  be  placed  farthest  from  the  target  (to 
convert  the  highly-monodirectional  incident  beam  into  a  more  dispersed  beam) 
and  subsequent  layers  would  be  included  for  gamma-ray  (and  neutron)  absorp¬ 
tion,  should  be  effective.  The  dependence  on  lateral  dimensions  and  non- 
idealized  geometries  has  not  yet  been  investigated. 

The  main  conclusions  to  be  drawn  are  the  following:  (1)  target  configura 
tion  plays  a  major  role  in  determining  shield  effectiveness;  (2)  apparently 
minor  geometric  effects  (e.g.,  transferring  small  thicknesses  from  one  layer 
to  another)  can  have  surprisingly  significant  effects  on  target  dose;  and 
(3)  the  dose  distribution  over  a  large  target  volume  is  not  uniformly  modified 
by  layering,  l.e.,  a  shield  mass  capable  of  achieving  only  modest  overall  dose 
reduction  may,  by  geometric  means,  provide  adequate  shielding  for  particular, 
vulnerable  locations.  The  objectives  of  the  Phase  I  study:  (1)  to  develop  a 
suitable  formalism  and  (2)  to  demonstrate  the  importance  of  geometric  optimi¬ 
zation  in  a  limited  number  of  moderately  simple  cases,  have  been  fully 
realized.  In  a  Phase  II  submittal  it  is  proposed  to  apply  the  principles 
already  demonstrated  to  much  more  complex  designs  which  can  be  used  in  mixed 
neutron-gamma  fields  and  in  cases  where  shielding  is  required  over  a  solid 
angle  approaching  a  half-space. 
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APPENDIX  A 


INTERSECTIONS  OF  GAMMA-RAY  PATHS  WITH  SHIELD  AND  TARGET  BOUNDARIES 
Ray  to  Plane 

The  planes  are  all  x  =  constant  surfaces  with  r  <  R,  where  R  Is  the 
radius  of  the  shield  layer  or  target  cylinder.  If  a  =  0,  the  ray  is  parallel 
to  the  plane  and  there  Is  no  Intersection.  Otherwise,  we  set 

xc  =  X-f 2  ,  a  >  0  (A-la) 

xc  =  Xn  ,  a  <  0  ,  (A-lb) 

where  i  is  the  number  of  the  shield  layer  that  the  photon  position  (x,y,z)  is 
in.  Then,  we  must  find  the  point  (xc,yc,zc)  where  the  line 

*c  '  x  Vc  “  y  zc  -  z 

- =  -  =  -  (A-2) 

a  B  7 

intersects  the  plane 

x  =  xc  ,  (A-3) 

which,  for  a  #  0,  is  obviously 


B 


yc  =  y  +  -  (*c  -  x) 

a 

(A-4a) 

7 

zc  *  z  +  -  (xc  -  X)  . 

(A-4b) 

a 


We  also  calculate  the  radial  distance  to  the  intersection  point, 

rc  =  tVc2  ♦  zc2]l/2  (A-5) 

and  check  to  see  if  rc  <  R.  If  not,  the  photon  intersects  the  cylindrical 
surface  before  reaching  the  planar  surface. 

Ray  to  Cylinder 

The  cylindrical  surfaces  are  of  the  form 


and  pick  the  solution  for  which  (yc  -  y)/B  >  0.  We  then  compute 


*c  =  *  +  “  (Yc  -  y) 

B 


(A-16a) 


zc  =  z  +  -  (yc  -  y) 

B 


(A-16b) 


If  r  =  0  and  B  =  0,  there  Is  no  solution.  If  xc  >  or  xc  <  X^,  the 
ray  exits  the  planar  ends  before  reaching  the  cylindrical  surface. 

Ray  to  Sphere 

A  sphere  of  radius  Rs  centered  at  ( Xs ,0,0)  can  be  represented  by 

(xc  -  h)2  +  Yc2  +  zc2  =  Rs2  •  (A‘17) 

For  all  cases  of  Interest  to  us,  x  will  be  less  than  xc,  since  the  shield 
layers  are  all  to  the  left  of  the  target  central  sphere.  If  «  s  0,  there  will 
thus  be  no  Intersection. 

If  a  >  0,  we  find  again  two  possible  solutions, 


-B  t  [62  -  4AC]V2 


( A-18 ) 


where 


A  =  — 
.2 


( A— 19) 


[By  +  yz  B4  +  r4 
- — 

ot  a4 


x  -  X< 


(A-20) 


(By  +  yz  \  B^  + 

-  J  +  -  x2  _  rs2  . 

a  /  <*2 


(A-21) 


If  Er  -  4AC  <  0,  there  are  no  solutions  and  the  ray  does  not  Intersect  the 

p 

sphere.  If  8  -  4AC  =  0,  there  Is  only  one  solution  and  the  ray  Is  tangent  to 

O 

the  sphere.  If  B  -  4AC  >  0,  there  are  two  solutions  and  we  take  the  one  for 
which  xc  -  Is  smaller.  Then  we  find 


yc  =  y  +  “  (*c  -  x) 


(A-22a) 


ZC  =  Z  +  -  (*c  -  x) 


( A— 22b ) 
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APPENDIX  B 

AVERAGE  ENERGY  LOST  IN  COMPTON  COLLISION 

There  is  a  one-to-one  correspondence  between  the  energy  of  the  scattered 
photon,  hv',  and  the  scatter  angle,  0,  in  the  Compton  collision  of  a  photon  of 
energy  hv  and  a  free  electron;  according  to  Evans  [16J  this  relation  can  be 
written 


hv'  = 


(B-l) 


1  -  cose  + 


where  nuc^  =  0-511  MeV  is  the  electron  rest  mass  energy.  We  let  n  be  the 


incoming  photon  energy  in  units  of  m0c  » 


(B-2) 


and  define  the  reduced  energy  variable 


Then  Eq.  B-l  can  be  written  in  the  form 
6=1+  n(l  -  cose) 


or  its  inverse 


n  +  1  -  « 


cose  = 


(B— 3) 


(B-4) 


(B-5) 


The  Klein-Nishina  differential  scattering  formula  for  unpolarized 
incident  radiation  is  given  by  Evans  [16]  as 

rr>2  /V'  v2  T  v  v'  o  1 


d(0a) 


=  !£  ^l2  f---- 


2*  sine  de  . 


(B-61 


where  r0  is  the  classical  electron  radius.  Equation  B-6  gives  the  probability 
that  a  photon  of  energy  hv  will  scatter  within  de  about  0  with  final  energy 
hv'  and  has  the  form  d{eo)  =  f(e)  sine  de.  We  can  rewrite  Eq.  B-6  in  terms  of 
6,  obtaining 


d(e*)  = 


"o2  r  i 

Hfi  +  r 


sin  e  I  sine  de  . 


(B-7) 


Because  of  the  correspondence  between  6  and  e,  given  n*  expressed  in  Eqs.  B-4 
and  B-5,  we  can  write 


d(eor)  =  f(e)  sine  de  =  9(4 1 r»)d6  , 

from  which  we  obtain 

de  dcose 

g(6|n)  =  f ( © )  sine  —  =  -f (e )  - 


(B-8) 


(B-9a) 


Equation  B-9b  is  a  restatement  of  the  Klein-Nishina  scattering  formula  which, 
when  normalized,  gives  the  probability  density  function  for  s  given  n. 


The  energy  lost  In  a  Compton  collision  Is  transferred  to  the  struck 
electron;  the  energy  of  the  secondary  electron  is  thus  [16] 


T  =  hv  -  hv'  , 

which  can  be  written  in  terms  of  6  as 


(B-lOa) 


T  =  hv 


H> 


(B-10b) 


We  want  the  average  energy  lost  per  Compton  collision,  which  is  given  by 


/  T(6)g(8|n)dn 

«0 _ _  N 

0 

/  g(«ln)dn 


(B-ll) 


.  4-f ( 


l\> 


where,  by  Eq.  B-4 


«0  =  «(8  =  0)  =  1  ( 
and 

6 1  =  S  { e  *  ir )  =  1  +  2n  ( 

Performing  the  integrations  in  Eq.  B-ll,  we  obtain 

hv  ir02  u2  _  2n  -  3  4  2n2  -  2n  -  2  2n(l  +  n) 

N  =  - < - s -  jt-n ( 1  +  2n)  + - +  - 5 

n  (  n  n  n(l  +  2n)  (1  +  2n ) ^ 

2n(4n2  +  6n  +  3)1 

3(1  +  2n)3  j  ( 

and 

Ttr02  (n2  -  2n  -  2  4  2n(l  +  n)) 

D  =  - < - * - J.n(l  +  2n)  +  -  +  - >  .  ( 


APPENDIX  C 


METHOD  TO  FORCE  SCATTER  TOWARD  CENTRAL  SPHERE 


A  photon  with  Incoming  direction  8  -  (a,a,y)  interacts  in  the  shield  at 
position  (x,y,z).  We  wish  to  force  a  scatter  Into  the  solid  angle  subtended 
by  the  target  sphere,  centered  at  (Xs,0,0)  and  having  radius  Rs,  as  shown 
in  Fig.  C-l.  We  first  find  the  direction  Q0  =  (a0>Bo»*o)  ^rom  (x>y»z)  t0 
(Xs,0,0),  which  is  given  by 


Xs  -  x 


a0 

(C-la) 

h 

-y 

Bo  =  “ 

(C-lb) 

h 

-z 

r0  =  —  > 

( C-lc ) 

h 


where 


h  =  [(Xs  -  x)2  +  y2  +  z2]1/2  .  (C-2) 

We  then  find  the  scatter  angle,  e0,  between  8  and  80,  which  Is  given  by 
ra(Xs  "  x)  ~  By  “  *z) 


en  =  cos 


r 


(C— 3) 


since  cose0  =  a  ■  80. 

The  sphere  projects  a  circle  of  radius  Rs.  In  order  to  force  a  scatter 
into  the  solid  angle  of  the  sphere,  we  must  pick  a  scatter  angle,  e,  in  the 
Interval  ( e0_ , e0+ ) ,  where 

®ot  =  eo  1  9m  >  (C-4) 

with 

-1 

em  *  tan  —  .  (C-5) 

h 
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}  >  r 

V-J 


This  Interval  defines  the  region  between  two  cones  that  are  tangent  to  the 
sphere.  To  force  e  Into  this  region,  we  pick  cose  uniformly  from  that  region 
for  isotropic  scattering  or  from  the  Klein-Nishina  formula  for  Compton 
scattering.  Thus  for  isotropic  scattering,  we  pick 


cose  =  cose0_  +  £[cose0+  -  cose0_] 


^ith  weight 


[cose0+  -  cose0_] 


>r  Compton  scattering,  we  pick  s  from  g(6|n)  with 
«0  =  1  +  n[l  -  cose0_] 

=  1  +  n[l  -  cose0+] 
a  id  apply  the  weight  factor 

«1 

/  g( 5 | n)dn 
«0 


where  n,  g,  and  0  are  given  in  Appendix  A.  Then  we  obtain  e  from 


( C— 6 ) 


(C-7) 


(C-8a) 

(C-8b) 


( C— 9 ) 


n  +  1  -  s' 


e  =  cos 


(C-10) 


Now,  given  e  we  need  to  force  $  into  an  interval  2A<J>  wide,  where 


_!  f(M  -  )1/2 ‘ 


A*  =  tan 


u  =  h  tan(e0  -  e) 

The  weight  factor  for  forcing  $  into  an  interval  of  width  2a*  Is 


(C-ll) 


(C-12) 


APPENDIX  D 


MCNS  CODE  INPUT/OUTPUT 

The  results  reported  herein  were  generated  using  a  FORTRAN  code  called 
MCNS,  whose  logic  is  outlined  in  Chapter  III.  Input  to  the  code  is  fairly 
simple;  the  required  input  format  is  identified  in  Table  D-l  and  the  input 
variables  are  defined  in  Table  0-2.  The  code  uses  several  logical  variables 
to  specify  available  options.  For  instance,  the  value  LISO  =  True  causes  all 
scattering  events  to  be  treated  as  isotropic  and  nonenergy-degrading,  which  is 
useful  in  code  verification;  the  normal  value  LISO  =  False  allows  the  code  to 
sample  post-Compton  conditions  from  the  Klein-Nishina  differential  scattering 
formula.  Also,  the  logical  variable  LEVEN  is  used  to  specify  whether  or  not 
the  simplified  geometry  of  Chapter  II  is  being  used.  If  it  is,  the  code  input 
requirements  are  simplified. 

The  code  operates  in  double  precision  and  at  present  is  dimensioned  to 
allow  20  shield  layers,  8  energy  Intervals  in  the  Input  spectrum,  10  annular 
rings  for  scoring  the  target  dose  and  3  elements  (hydrogen,  oxygen  and  lead). 
The  code  Is  reasonably  efficient,  due  to  the  use  of  last-flight  estimators, 
analytic  equivalence  and  the  forcing  of  scattered  photons  into  the  solid  angle 
of  the  target  sphere.  Typically,  8,000  incident  photon  histories  were  used 
for  a  shield  of  four  mean  free  paths  thickness  and  2,000  histories  for  shields 
whose  total  thickness  was  one  mean  free  path.  Sample  standard  deviations  were 
typically  less  than  0.5  percent. 

The  code  was  written  in  modular  form  and  thus  can  easily  be  extended  to 
consider  more  complex  situations.  For  instance,  It  would  be  a  straightforward 
matter  to  replace  the  plane-parallel  layers  by  conic  or  curved-surface  layers, 
to  add  new  elements  to  the  cross-section  'ibrary  or  to  allow  layers  whose 
densities  or  compositions  varied  from  on.’  to  another. 


61 


TABLE  0-1.  MCNS  CODE  INPUT 


Record 

Number  Variables  Format 


1  I DATE ( I ) ,  I  =  1,3  3110 

2  IW,  IX,  IY,  IZ  4110 

3  LRITE,  LRAL,  LISO,  LEVEN,  LOUT  5(9X,L1) 

4  NHIS,  NSH,  NUME,  NE,  NRING  5110 

I  NHIS  <  0,  Stop 

5  RSAT,  TSAT,  RSPH,  RCEN,  CUTOFF  4F10.0 

If  LEVEN  =  True,  go  to  record  9 

6  SPSH(I),  I  =  1,2,...,  NSH  7F10.0 

7  TSH(I),  I  =  1,2,...,  NSH  7F10.0 

8  XSAT  F10.0 

Go  to  10 

9  CAPT,  CAPD,  SMD  3F10.0 

10  WTS( I ,K) ,  K  =  1,2,...,  NE;  I  =  1  7F10.0 

Repeat  record  10  for  I  =  1,2,..,,  NSH+1 

11  RHO(I),  I  =  1 ,2, . . . , NSH+1 
If  NUME  >  0,  go  to  record  13 

12  ES  F10.0 

Go  to  record  4 

13  EC ( J > ,  J  =  1,2, .. . ,NUME+1  7F10.0 

14  SE( J) ,  J  =  1,2,. ...NUME  7F10.0 

Go  to  record  4 


TABLE  0-2.  MCNS  CODE  INPUT  VARIABLE  DESCRIPTIONS 


Variable 

Name 


Description 


IDATE(I), 
1=1, 2, 3 

IW,IX,IY,IZ 

LRITE 


LEVEN 


NRING 


CUTOFF 


Month,  day,  run  number 

Random  number  generator  initial  seeds 

Logical  variable;  if  T,  detailed  results  are 
output;  in  F,  only  final  results  are  output 

Logical  variable;  if  T,  Rayleigh  scattering 
is  simulated;  if  F,  Rayleigh  scattering  is 
neglected 

Logical  variable;  if  T,  scattering  is 
isotropic  and  no  energy  is  lost  in  a  scatter; 
if  F,  scattering  follows  Klein-Nishina  model 

Logical  variable;  if  T,  shield  layering  is 
even  and  simplified  input  option  is  invoked; 
if  F,  detailed  input  is  used 

Logical  variable;  if  T,  results  also  written 
to  output  file  on  logical  unit  4;  if  F, 
results  only  written  to  line  printer 

Number  of  MC  histories 

Number  of  shield  layers 

Number  of  energy  intervals  in  source 
spectrum 

Total  number  of  elements  in  all  shields  and 
the  target 

Number  of  concentric  rings  on  first  shield 
layer  for  scoring  contributions  to  dose 

Radius  of  target  cylinder 

Thickness  of  target  cylinder 

Radius  of  target  sphere 

Radius  of  central  sphere  for  forced  target 
scoring 

Cutoff  value  of  the  product  of  photon  energy 
and  history  weight  below  which  the  photon 
history  is  terminated 


TABLE  D-2.  MCNS  CODE  INPUT  VARIABLE  DESCRIPTIONS  (Continued) 


Variable 

Name 


Description 


Units 


SPSH(I) 

TSH(I) 

XSAT 

CAPT 

CAPD 

SMD 

WTS(I.K) 

RHO(I) 

ES 

EC(J) 

SE(J) 


Spacing  between  shield  layers  1-1  and  I 

(cm) 

Thickness  of  shield  layer  I 

(cm) 

X-coordinate  of  target  left  face 

(cm) 

Total  thickness  of  all  shield  layers 

(cm) 

Distance  from  the  front  face  of  first  shield 
layer  to  front  face  of  target  cylinder 

(cm) 

Separation  between  front  face  of  target 
cylinder  and  back  face  of  nearest  shield 

1  ayer 

(cm) 

Weight  fraction  of  element  K  in  shield 
layer  I  >  0,  or  target  (1=0) 

Density  of  shield  layer  (I  <  NSH)  or 
target  (I  =  NSH-t-1) 

(gm/cm3) 

Source  energy,  if  monoenergetic  (NAME=0) 

(MeV) 

Energy  spectrum  cutoff  points  such  that 

Jth  energy  interval  is  between  EC(J-l) 
and  EC ( J ) 

(MeV) 

Fractional  part  of  spectrum  in  energy 
interval  J 
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